Spin transfer torque in continuous textures: 
semiclassical Boltzmann approach 
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We consider a microscopic model of itinerant electrons coupled via ferromagnetic exchange to 
a local magnetization whose direction vector n{r,t) varies in space and time. We assume that 
to first order in the spatial gradients and time derivative of n{r,t) the magnetization distribution 
function /(p, r, t) of itinerant electrons with momentum p at position r and time t has the Ansatz 
form f{p,r,t) = /||(p)n(r,t) + /ir(p)n x Vrn + /2r(p)Vrn + /it(p)n x dm + f2t{p)dtn. Us- 
ing the Landau-Sillin equations of motion approach we derive explicit forms for the components 
/|| (p)i /ir(p)) /2r(p), /it(p) and /2t(p) in "equilibrium" and in out-of- equilibrium situations for (i) 
no scattering by impurities, (ii) spin-conserving scattering and (iii) spin non-conserving scattering. 
The back action on the localized electron magnetization from the out-of-equilibrium part of the two 
components /ir,,f2r constitutes the two spin transfer torque terms. 
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PACS numbers: 72.25.Ba, 72.10.-d 



I. INTRODUCTION 



Recent experiments in spin-valve nanopillars^^, point 
contacts? and ferromagnetic nanowires^i^iii^ have 
demonstrated the possibility to "manipulate" the mag- 
netization by applying an electrical current instead of 
an external magnetic field. It is believed that this phe- 
nomenon might give rise to many technological appli- 
cations, (MRAM, fast magnetic switching, high quality 
microwave sources), provided that the current density 
necessary to manipulate the magnetization can be dras- 
tically reduced. 

This possibility of electrical current induced magneti- 
zation manipulation was predicted already ten years ago 
by J. C. Slonczewski and L. BergerSjlS. 

Nowadays, the consensus is that, in presence of an 
electrical current, the standard Landau-Lifshitz-Gilbert 
(LLG) equation that describes the dynamics of the mag- 
netization is modified. For ferromagnetic wires in which 
the magnetization texture M(r, t) — Ms'n{r, t) has a time 
and space dependent direction n(r,t), the modified LLG 
equation reads-W: 

dtu = 7Bcff (n) X n + an x dtfi — u • Vrn + (3n x u • Vrn. 

(1) 



where Boff(n) 



_i 5£{n) 



is the effective magnetic field 



Sn 

derived from the magnetic free energy and u, (3 are two 
phenomenological parameters where u is proportional to 
the electrical current density and polarization, and a is 
the usual phenomenological Gilbert damping parameter. 
The contribution u • Vrn is called the adiabatic term 
and can be derived from an additional term in the mag- 
netic free energy that takes into account the coupling 
to the electrical currenti^. The /3 term is in contrast 
non adiabatic and appears to play a role similar to the 
Gilbert damping term (and one conception of the spin- 
transfer torque introduction in the LLG equation even 
predicts (3 = a exactly^i^) . Recent micromagnetic numer- 
ical experiment o-'^^d^d^d'^ with this modified LLG equa- 



tion have greatly clarified the qualitative roles played by 
the two phenomenological parameters u and (3. In partic- 
ular, it has been proved that in the absence of the (3 term 
there is no current induced steady domain wall motion 
below a finite (very high) critical current density Uc. For 
a non zero j3 term, in the simple case of a perfect sam- 
ple it can be shown that the speed of the domain wall 
is {f3/a)u. More quantitatively, results compatible with 
the various experiments are obtained for a ratio (3/a > 1 
equivalent to f3 ~ 10~^. 

The derivation of the LLG equation above usually 
rests on a two steps argument. It is assumed that an 
itinerant ferromagnet can be modeled as a two "com- 
ponents" system: (i) non moving and ferromagnetically 
ordered electrons (called hereafter d) described by a clas- 
sical magnetization vector n(r, t) that varies slowly with 
time t and position r; (ii) current-carrying itinerant elec- 
trons (called s) coupled to the d electrons via a ferro- 
magnetic exchange energy Asd > such that the ef- 
fective one-electron quantum Hamiltonian has the form 



.y-^/ _|_ 



n[r, t) ■ a , where / is the 2x2 
is the vector of Pauli matrices 



identity matrix and a 
with eigenvalues ±1. 

The first step consists in finding the quantum 
average itinerant electron magnetization m(r, t) — 
— /iBTr{/3(n) o- } as a function of the " quasistatic" 
n(r, t) {p{n) is the itinerant electron density matrix that 
depends on n{r,t)). In a second step one substitutes 
the resulting m(r,t) back into the LLG equation of the 
localized electrons magnetization: 



dtTi = 7Bcff(n) X n -I- an x dtn - 



1 



MsUB 



Asdm{r,t) X n. 

(2) 

For domain walls, Zhang and Li (ZL)^^ were the first 
to make a transparent derivation along this line of rea- 
soning. They found that in out-of-equilibrium situation 
the above back action not only produces the two spin 
transfer torque terms but also leads to correction of the 
Gilbert damping term and gyromagnetic ratio. Prior to 
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ZL, Zhang, Levy and Fert— had obtained corresponding 
resuhs for spin- valves. Both works however rest on a phe- 
nomenological equation of motion for the magnetization 
m{r,t) of itinerant electrons where key ingredients are 
put by hand, especially spin flip scattering time and adia- 
baticity of the itinerant electron spin current with respect 
to local magnetization n{r,t). Beside these phenomeno- 
logical descriptions, various microscopic approaches have 
also flourished in the last years. For spin-valves sys- 
tems several groups used a scattering matrix approach 
(for a review see^^). For ferromagnetic wires the con- 
cept of local spin reference frame was often invoked so 
as to exhibit a direct coupling between spin current of 
itinerant electrons and gradient of the local magnetiza- 
tion n(r, To the best of our understand- 
ing however, none of these different works really succeeds 
to establish the modified LLG equation as written above. 

Very recently two independent works, by Tserkovnyak 
et al^ and Kohno et al^^ , based on different theoretical 
techniques presented a direct microscopic derivation of 
the two additional spin torque terms. Both works show 
that the /3 term requires the existence of a spin flip like 
scattering mechanism (e.g. spin non-conserving scatter- 
ing like spin orbit, magnetic impurities ...). They show 
that this mechanism is also responsible for the appear- 
ance of an effective a Gilbert damping term induced by 
the itinerant electrons. In fact Tserkovnyak et at argue 
further that one should find a = f3 for itinerant ferro- 
magnetic systems where the magnetism comes from the 
exchange interaction between the itinerant electrons (e.g. 
the effective local n(r, t) is in fact the itinerant electron 
magnetization m(r, t) itself). The results of Kohno et al. 
were obtained from diagrammatic linear response theory 
and concerned only "integrated" physical quantities like 
the local magnetization m(r, t). In contrast Tserkovnyak 
et al. used the Keldysh quasiclassical Green function 
technique that in principle could allow determining the 
full magnetization distribution function /(p, r, t), and 
thus might give a deeper understanding of the system. 

In this paper we reconsider the entire problem within 
a Boltzmann approach. This provides an intuitive and 
hopefully pedagogical semiclassical picture of the equa- 
tions of motion of the charge and spin distribution func- 
tions of the itinerant electrons, in a space-time dependent 
magnetization field n(r, t). Our main assumption is that 
around any time space position the direction of n(r,t) 
can be arbitrary but its gradients VrW and dtn must be 
slow enough so that only terms parametrically linear in 
these gradients are important (e.g terms like VrW, dtn, 
n X Vpn and n x dtn). This parametrization is different 
from that of Kohno et al and Tserkovnyak et al. which 
is in fact a linear theory around the uniform magnetiza- 
tion case (e.g. n(r,t) — z-|-u(r,t) with |u(r,t)| ^ 1). 
The difference might appear subtle, but it leads to dis- 
tinct properties already in the equilibrium situation as 
compared to Tserkovnyak et al.. 

The Boltzmann method provides a complementary 
physical picture that is intermediate between the purely 



phenomenological macroscopic equations of motion ap- 
proach of ZL on one side and the less intuitive but more 
microscopic quantum linear response or Keldysh meth- 
ods on the other side. More precisely, on the one hand, 
since we deal directly with the charge and spin distribu- 
tion functions, we have a deeper understanding than that 
provided by ZL method. We have also access to more mi- 
croscopic physical quantities, and can in principle study 
general AC and thermoelectric effects not accessible to 
ZL method^i. On the other hand, the physical picture 
provided by the Boltzmann method is more transparent 
than the microscopic approach, it also allows understand- 
ing more clearly at which level spin-flip scattering differs 
qualitatively from spin- conserving scattering. Lastly, the 
Boltzmann method provides the natural framework to 
see where quantum corrections could be important. In- 
deed, in a forthcoming paper we will reexamine our re- 
sults using the Keldysh Green function method in the 
quasiclassical approximation. This constitutes the natu- 
ral theoretical framework to build from first principles the 
equations of motion for the distribution function in the 
presence of elastic and inelastic collisions. The Keldysh 
approach appears necessary because the construction of 
the collision integral in the Boltzmann picture is purely 
phenomenological so that, on a more microscopic level, it 
is not clear if there are important quantum and gradient 
corrections to the Boltzmann collision integrals. 

The paper is organized as follows. In the first section, 
going back to textbook a^^i^^ we rederive the collisionless 
transport equations for the charge and magnetization dis- 
tribution to first order in time-space gradients. We show 
that without any collision the magnetization distribution 
function /(p, r, i) is not coUinear to n(r,t) and that, in 
particular, there is a contribution coUinear to n x 
that gives rise to a finite equilibrium spin current. In the 
next two sections with study the influence of elastic scat- 
tering by impurities. Assuming Boltzmann type collision 
integrals, we first consider the effect of spin-conserving 
collisions, in equilibrium and in out-of-equilibrium situ- 
ation. Already at equilibrium we obtain surprising re- 
sults for the components of the distribution that are not 
coUinear to n(r, t). We show in particular that the equi- 
librium spin current is rotated in the direction Vr'T- by 
an angle that depends on the ratio between the elastic 
scattering time and the effective Larmor time. In out- 
of-equilibrium situation we find correspondingly a com- 
ponent of the magnetic distribution function coUinear to 
Vrn that could in principle give rise to a /3 term as back 
action after p momentum integration. (Un)fortunately 
it vanishes after p momentum integration as spin conser- 
vation imposes. We then consider collisions that lead to 
spin flip. For the case of a uniform magnet, starting from 
the known form of the collision integral of each eigen-spin 
distribution, we build a collision integral invariant under 
spin basis change. Extending phenomenologically this 
collision integral to the case of a non uniform ferromagnet 
we describe how the equilibrium and out-of-equilibrium 
properties are modified by the spin Uip scattering. In 
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particular, from the different components of the itinerant 
electron magnetization we evaluate their back action on 
the d electron magnetization and then extract explicit ex- 
pressions for the induced Gilbert damping a2t, modified 
adiabatic torque term u and finally the f3 torque term. 
Our expressions for a2t and P coincide exactly with ZL 
approach and to leading order also coincide with results 
of Kohno et al for spin- isotropic spin flip scattering, when 
appropriate changes of notations are made. Concerning 
the parameter u, that is purely phenomenological in the 
ZL approach, our result is similar to that of Kohno et al. 
Note finally that, in order to improve reading, most of the 
calculations concerning spin flip scattering are put in the 
Appendix whereas the main text essentially provides the 
physical picture emerging from the modified distribution 
functions. In a last concluding section we discuss pos- 
sible extensions of our approach to itinerant ferromag- 
nets or ferromagnetic Fermi liquids^^'^^ and spin-valve 
systems^. 



II. SEMICLASSICAL TRANSPORT THEORY 

A. Model and Ansatz 

The effective one-electron quantum Hamiltonian of the 
itinerant free electrons coupled to the localized electron 
magnetization is of the form 

Hir) = [-^Vl + Vir)]i + ^n{r,t)- a, (3) 
2m 2 

where V{r) = — eEr is the potential induced by an ex- 
ternal uniform electric field. The intrinsic difficulty to 
find the equilibrium and out-of-equilibrium density ma- 
trix p{r,r',t) associated to this Hamiltonian originates 
from the non commutation of the "Zeeman" term with 
the kinetic term, due to the spatial variation of n(r,<). 

However, as we consider here domain walls where the 
characteristic length of the magnetization gradient is 
large (10-100 nm) compared to the electron mean free 
path, quantum transport is not pertinent and electron 
diffusion is a more appropriate framework. Therefore we 
simplify the problem by assuming that the spatial degrees 
of freedom r, p are classical commuting variables and not 
quantum operators, and retain only the non trivial com- 
mutation rules of spin degrees of freedom. The effec- 
tive semiclassical Hamiltonian of the itinerant electrons 
is thus 



H(r,p) = [i^ + y(r)]/- 
2m 



sd 



n(r,t) 



(4) 



where we have denoted by usgii — the effective Larmor 
frequency. We further define Tsd = and tsd — vpTsd 
the Larmor time and Larmor length respectively {vp is 
the Fermi velocity). 



The semiclassical quantity that corresponds to the den- 
sity matrix is now the spinor distribution function 



1 



/(p,r,i) = -[/(p,r,i)/ + /(p,r,<)- a 



(5) 



The physical quantities such as local particle density 
n(r,f), particle current density j{r,t), magnetization 
density m(r, t) and spin-current tensor density J(r, t) 
are obtained by integration on these distributions: 



where 



n(r, t) 

m(r, t) 
J(r,i) 



dr : 



/dTTr,{/(p,r,t)} 
/dT/(p,r,t), 
/dr^/(p,r,t)^ 
-^B/dTTr<,{/(p,r,t) & } 
-fiB J dTf{p,r,t) 
-Ms /dT£-/(p,r, t) 



(6) 



dp 



with p the unit vector of direction p, ep = ^ and i^(e) 



the 3D free electrons density of states. 



\/2m^e 

Our main assumption is that around any space-time 
position r,f, the direction n{r,t) can be arbitrary but 
its gradients Vrn (resp. dtfi) must be slow enough com- 
pared to the Larmor length i^d (resp. Larmor time) so 
that only terms linear in these gradients are important. 
Linearization in these gradients implies that time and 
space dependencies of the matrix distribution function 
/(p, r, t) are expanded on the possible directions n{r,t), 
n X Vrfi, n X dtTi, VpTi, and dtn. The Ansatz form we 
assume for /(p,r,t) compatible with this approximation 



/(p,r,i)= /(p), 
/(p,r,i) = /||(p)n(r,i) 

+(-sd [/lr(p)«. X VrTl + /2r(p)Vrn] 

+Tsd [/it(p)w X dtn + /2t(p)9tn] . 



(7) 



Note that (n, dtn, n x dtn) and (n, VrW, n x Vrn) con- 
stitute two distinct orthogonal bases for any spin vector. 
Thus a priori our Ansatz contains some redundancy since 
dtn, n X dtn are linear combinations of Vrn, n x V^n 
and reciprocally. The main reason why our Ansatz is 
nevertheless appropriate is that, in this extended ba- 
sis {n,dtn,n x dtn,Wrn,n x Vrn), to leading order, 
each component of the spin distribution is stationary (for 
DC field), space independent, and depends only on p"^-. 
Had we chosen the basis (n, dtn, n x dtn) to expand the 
spin distribution, each component would still be position- 
dependent to leading order. 

Note that with our normalization the components /ir, 
/2r, fit and /2t do have the same physical dimension 
as /||. As a consequence, the quantity mir(r, i) = 
—^iB J drfirij), r, t) has the units of a magnetization den- 
sity. 
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Following Landau-Sillin^^, to first order in the gradi- 
ents of the distribution function, the " Liouville" equation 
of motion of the distribution is obtained from 



dl 
dt 



1 



\{f.H}-\{Hj}=i[h (8) 



where [A, i?]^ denotes the commutator (resp. anticom- 
mutator) of ^ and B and {A, B} = VrAVpB- VpAVr^ 
designates the classical Poisson brackets. represents 
the spin operator for the collision term. In the absence 
of collision, we obtain the following coupled set of equa- 
tions of motion for the particle and spin components of 
the distribution function: 



(9t + ^Vr + eEVp)/ -^Vrn-Vp/ = 0, 
(5t + ^Vr + eEVp)/ -%iVrnVp/ 

-LOsdn X / = 0. 



(9) 



In these equations the symbols x and • mean vector prod- 
uct and scalar product of spin-vectors. For spatial quan- 
tities written in bold an implicit scalar product is un- 
derstood, namely ^Vr = Y,i m^''; ' ^p-^ = 
X^i ^rtn-Vp-f . For later use, we further write / = f°+g 
to separate the "equilibrium" {E — 0, but n may depend 
on t) contribution /° from the out-of-cquilibrium contri- 
bution g; and accordingly for each component. 



B. equilibrium properties 

In the case of no applied electric field, the above equa- 
tions have simple solutions. 

a. Homogenous ferromagnet: n{r, t) = no In the 
absence of an electric field and for a time and space in- 
dependent direction n(r,t) = no the stationary distribu- 
tions are well known: 



/0(p) = i(n^(e+)+n^(e-)), 
/(p) - /|°(p)no, 



with 

/|?(P) = hMe+) 



(10) 



"F(ep)), 



2m 



and nF(e) = (e'^^'"^) -f 1)-^ is the 



where 

Fermi statistic. The physical quantities are then: 

n(r,i) = ne = /dei.(e)/0(e), 
m(r, t) = m||no, 
with 

m|| = [iBPn^ = -hb J dei^(e)/|^(e). 



(11) 



that defines the polarization P. Note that the scalar 
functions /" and f^^ depend in fact of the energy e. 

b. Inhomogenous ferromagnet n(r, t) The naive ex- 
tension to a space and time dependent n(r, t) would be an 
unchanged particle distribution and a magnetization dis- 
tribution that follows n(r, i) adiabatically: f{p,r,t) — 
/|'j'(p)n(r, t). (Un)fortunately it is not the solution of 



the equations of motion above. Instead, to linear order 
in the space-time gradients, the solution is of the form: 

/(p,r,t) = /|^(p)n(r,t) 

+£sdfirip)n X Vrn + Tsdfitip)n X dtn 

(12) 



with 



and 



f?tiP) = -/f (P), 



.lai P. f , (f) 



/^(e)-/f(e)-^aj°(e) 



(13) 



(14) 



The above defined function (e) will appear many times 
in this paper and constitutes a key ingredient of most 
transverses magnetic quantities. The figure below plots 
this function /_L(e) compared to /^(e). 




FIG. 1: The two functions /_L(e) and f\\{^) for parameters 
ksT = 0.025meV and huJsd = Q.2eV . 

The physical consequences of these non adiabatic com- 
ponents are the following. The itinerant electrons mag- 
netization is m(r,t) = m||n(r,t) + Tsdmit'n{v,t) x dtn 
with rriit = As already pointed out by Zhang 

and Li^^, in the s — d picture, the back action from 
the component mu renormalizes the 7 term of the d 
electrons magnetization LLG equation. The component 
firiv) does not contribute to the magnetization because 
of its rotational symmetry in momentum space. However 
it contributes to a local equilibrium spin current tensor 
J{r,t) = £sdJirn{r,t) x Vrn with 



Ji) =-/.B/dr|^/0.(p) 



(15) 



dp. 



where we used the identity / -^p'p-' 



In the s — d 

picture, this equilibrium spin current has no back action. 
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In the itinerant ferromagnet picture where conduction 
electrons are responsible for the local magnetization it- 
self (e.g. Asd is the exchange interaction between the 
electrons), the divergence of this equilibrium spin cur- 
rent should be equivalent to the local exchange field term, 
as pointed out by several authors^i*^. This equilibrium 
spin current cannot be derived from the ZL macroscopic 
approach. As we will show in a forthcoming paper, it is 
also obtained with the quasiclassical Keldysh approach. 



III. SPIN-CONSERVING IMPURITY 
SCATTERING 

A. Collision integral and transport equations 

We now consider the scattering by impurities in a 
purely phenomenological manner. For spin-conserving 
scattering, a first guess is a local collision integral of the 
form 



J[/(p,r,t)] =li + X- & 



where 



= /rfT'wp,,p/(p',r,i)(l-/(p,r,t)) 
-Wp.p'/(p,r,t)(l - /(p',r,i)), 

(16) 

is the probability to scatter from momen- 
tum p to momentum p'. In the expression above 
(1 — /(p, r, i)) means (1 — /)/ —fa. In general 
/(p',r,t) can be non coUinear to /(p,r, i). As a conse- 
quence, due to the Pauli matrices properties, the above 
form leads to an unphysical imaginary contribution like 



«(wp',p Wp^p')(/(p',r,i) X /(p,r,t)) • <t . To prevent 
such a contribution it is necessary to antisymmetrize cor- 
rectly the previous expression: 



As a consequence, the magnetization in the presence of 
spin-conserving scattering and DC electric field driving 
the system out-of-equilibrium cannot have finite compo- 
nents for either m2rVrn or m2tdtn. Thus, in that situa- 
tion, as back action there are no induced a and f3 torque 
terms. We stress that this somewhat disappointing result 
concerning the nonexistence of a /3 term (finite m2r) for 
macroscopic quantities is not characteristic of what hap- 
pens to the microscopic distribution function. In fact, 
as we explain below, at the level of the spin distribution 
distribution function, in the presence of spin-conserving 
impurities in out-of-equilibrium situation, there is a finite 
energy resolved magnetization m2r(e) that results after 
angular momentum integration (or a finite angular mo- 
mentum resolved magnetic density tensor TO2r(p) that 
results after energy integration). So, the itinerant elec- 
tron magnetization is not far from having a finite m2r- 

Spin-conserving scattering by impurities produces 
other interesting results that we now describe in both 
equilbrium and out-of-equilibrium situations. With our 
general Ansatz linear in the gradients, we obtain the 
following set of coupled equations for each "stationary" 
component: 

eEVp/ = X, 
eEVp/[i =X||, 

eEVp/2r + jf^(^/|| + ^Vp/) + LJsdflr = l2r, 
eEVp/ir — tOsdf2r — llr, 

eEVp/2t + ^/|| + uosdht - 
eEVp/if - ujsdf2t = lit, 

(20) 

The components /, /y are completely decoupled from the 
others. By contrast, /ir,/2r are coupled and depend on 
/, /||, and /it, /2t are coupled and depend only on /|[. 



/(p',r,t),(l-/(p,r,<)) 
/(p,r,t),(l-/(p',r,<)) 



(17) 



With this form, the equations of motion for the charge 
and spin components of the distribution function become: 



(5t + ^Vr-FeEVp)/ 
(ft-f ^Vr + eEVp)/ 



^VrnVp/ =1, 
2^Vrn- Vp/ 

-UsdU X / = X, 



(18) 



B. equilibrium properties 

To go further we need to specify the form of the proba- 
bility Wp p' . To simplify the calculation we only consider 



the case of isotropic scattering Wr 



As usual we define the inverse scattering time t by: 



dp^ 

(27r7i)3 



{ep)S(ep-£p,) =w{ep)iy{ep) (21) 



with 1 — j (ir'?iip^p/(/(p') — /(p)) (similarly for X) and 
where we use the property that for elastic scattering 

Wp.p' = Ifp'.p. 

Whatever the form of Wp^pi , by integrating these equa- 
tions of motion over p the resulting continuity equations 
for the particle density n(r, t) and local magnetization 
m{Y,t) are the same as in the coUisionless situation in 
absence of electric field, namely: 



dtn{r, t) + divj(r, t) = 0, 
dtin{r, t) + divJ(r, t) — Luisdn^r, t) x m(r, t) = 0. 



(19) 



and ri{ep) = The ratio rj is anticipated to be much 
smaller than unity. Nevertheless, the calculations will 
be performed for a general rj. With this specific form 
and in the absence of an electric field the previous ex- 
pressions of /(p),/||(p),/it(p), and /2t(p) = are still 
solutions; their corresponding collision integral vanishes 
because all these distributions are function of ep only. By 
contrast, the previous form of /ir(p) is such that for a 
general Wp,p' (and in particular for isotropic scattering), 
the corresponding collision integral Xir does not vanish. 
It thus implies a non vanishing /2r(p) component that 
may modify /ir(p) as a back action. One finds that the 
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new self consistent equilibrium solutions fi^{p), /2r(p) 



are: 



/?r(P) 
/2°r(P) 



(22) 



where f± is as given in (jl4p but divided by 1 + ry^. 
The main modification induced by the spin-conserving 
scattering is thus that, in equilibrium, the spin cur- 
rent has now two components: a component coUinear 
to —n{r,t) X VrTi and a component coUinear to — Vrn, 
of smaller amplitude by a factor rj = The total 

modulus of the equilibrium spin current is smaller by 
a factor , ^ compared to the collisionless situation. 



Apart from this change of modulus, the scattering has 
thus induced a rotation in spin space of the equilibrium 
spin current vector. The ratio between the two trans- 
verse components of the current defines a rotation angle 
6 by ta,nd = i] = —. Let us stress that it is rather 
unusual to find any modification of the equilibrium dis- 
tribution by elastic scattering; it might be a hint that 
some "quantum" corrections have been neglected. 

From now on, to simplify the expressions, we assume 
that T(e) is constant and independent of energy. When 
necessary, we comment on the validity of our results for 
an energy dependent T(e). 



C. out-of-equilibrium properties 

Let us now consider the effect of the electric field. In 
the presence of an infinitesimal electric field, to each 
previous component fa{p) will be added an out-of- 
equilibrium component ga(p) = 5q(p, e)- To linear or- 
der in electric field and for isotropic scattering, standard 
calculations lead to the solutions 



g{p) = -eE^rajO(e), 
g||(p) = -eE^r9,/|^(6). 



(23) 



The components git,g2t are solutions of the following 
coupled equations: 

7^^511 + ^sd9it = 7 / ^(g2t(p', e) - 52t(p, e)), 



dp' 

eEVp/i° ~ujsdg2t = i/— (5l^(p',e)-.glt(p,e))• 



The solutions are immediate: 



.9it(p) = -5||(P), 

g2t(p) = 0. 



(24) 



(25) 



The last two components are solutions of the following 
coupled equations: 



= 7/^(.92r(p',e)-fe(p,e)), 

= 7 / ^(5ir(p': e) - 5ir(p, e))- 



The general solutions of these equations are linear in E 
and have the parametric form: 



with 



gij&n X 



(27) 



and similarly for .g2rVrM. For the quantities xi^2{£) and 
?/i,2(e) we obtain finally: 



X2 



l+ay-' 2e 



dji 



2rf 



yi = + 

2/2 = -'jxi 



n 

f± 



1 j-0 

r,J\l 



(28) 



Even if gir and g2r are not simple distributions, from 
their rotational symmetry properties in momentum space 
we deduce that there is no finite spin current associated. 
By contrast, it is not clear if any finite magnetization 
will survive after angular and energy integration. Thus, 
it appears interesting to define two partial magnetiza- 
tion densities for each components. On the one hand, 
the energy resolved transverse magnetization densities 
mir(e), m2r(e) are obtained after integration over angles 
alone, and on the other hand the angle resolved magneti- 
zation tensors TOir(p), "^2r(p) are obtained after energy 
integration. These quantities are defined through: 



mir(r, t) 



-fiB J dTgir{r,p,t) 
i^-^ J de miAe), 



= &4^/^-ir(P), 



dp 



(29) 



in 



and similarly for TO2r(r; t). Performing the angular inte- 
gration, we obtain the following expression for the energy 
resolved transverse magnetization densities: 



mir(e) 



W2r(e) = -rnj_(e), 
with 

m±{e) ^ -fiBi^{£){hdef± + f±), 

W||(e) ^ -flBl^{e)f^|■ 



m 



The form of TOx(e) makes it clear why J de TOj^(e) = 0. 
A possible very challenging experimental test of our ap- 
proach would be to measure, with a magnetic STM tip, 
these effective itinerant electrons energy resolved trans- 
verse magnetic densities as a function of energy. In fact 
the existence of a finite m2r(e) in the absence of any spin 
fiip scattering would be a proof of the non trivial compo- 
nent /j_ that is also at the origin of the equilibrium spin 
current. 

Performing only the energy integration gives rise to the 
two angle resolved magnetic tensors m2r(p) and mir(p): 



"^ir(P) = iTl$)^MS'' - 3p''p^") - ir 
m2r(p) = iq^mi(5''^"-3p^pJ"), 
with 

m_L = -/is / de v{e)f±{e), 
m|| = -muB J de i/(e)/||(e). 



(31) 
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to prevent contributions to w;pp'(i",t) that are linear in 
gradients of n{r,t). But with our phenomenological ap- 
proach there is no way to guess their specific form. This 
is in fact another key reason why a quantum approach 
using Keldysh Green function technique might be useful. 

In the collision integral shown above, due to the ap- 
propriate antisymetrization and for elastic scattering, all 
the terms that are quadratic in the distribution func- 
tions cancel. It then appears that in the equations of 
motion detailed in the Appendix, the effective spin flip 
collision integrals that appear for each component have 
^ the forms: = /rfr'i«+ ± <-,)(/(p') - /(p)) 

drgiiip) = / dr^.g||^(p). (32) and Xf± = /dr'i(«;^+ ± <-,)(/a(p') + Up)) for 

a=||,l,2. 

In the following we shall consider only isotropic scat- 
tering in p space {w^^,{e) = w^'(e)). Accordingly we 
deflnc: 



These last two expressions show why a further angular 
integration cancels all the contributions induced by the 
existence of the non trivial component f±{e)- 

All these results are, in fact, valid for any spatial di- 
mension provided we use the corresponding density of 
states and angular integration. They can also be gener- 
alized to an energy dependent scattering time with more 
complicated expressions for the distributions 51, 2r- From 
the spin continuity equation, we know that we should 
find mir = This is indeed the case because we 

have the equality 



P 



For an energy independent scattering time we further- 
more obtain: 



e\E\7 



Note that, although our results are non perturbative 
in 77, they are valid order by order in 77 taking care of 
the fact that the leading order is I/77. This is in contrast 
with the linear response results of Kohno et al. where 
the standard "leading order term" leads to an unphysi- 
cal magnetization component that is cancelled only 
when vertex corrections that constitute infinite order re- 
summation are carefully taken into account. 

We can partially generalize our results in two ways^: 
(i) define distinct intraband scattering probabilities Wp p, 
for energies e±. (ii) define anisotropic scattering proba- 
bility Wp^p'. 



IV. SPIN-FLIP IMPURITY SCATTERING 



1 



We further define: 



^f(e±)Kep±7j^.<i), 



(34) 



1 

^sf 






1 


1 \ 




f 


= - + — ' 






< 1, 



(35) 



A. Spin-basis invariant collision integral 



and correspondingly 77*^^, 77 with 77*^^ < 77. 



To extend our Boltzmann approach to spin-fiip scatter- 
ing, a key step is to find a spin basis invariant formulation 
of a collision integral that characterizes a spin-flip pro- 
cess. We explain in details in the Appendix how to find 
this collision integral in the case of uniform magnet, and 
how to phenomenologically generalize it to non uniform 
magnets. The resulting collision integral is: 



i«f = 1 / dr' 



<p'(r,i),[/(p',r,t),l-/(p,r,i)] + 
■7i^^p,(r,t),[/(p,r,i),l-/(p',r,t)] + 



with 

w: 



<"p')^ 

5«;'"<p')n(r,<) 



and 

/(p,r,t) = /(p,r,i)/-/||(p,r,t). & 

(33) 

The fact that /(p, r, t) appears in the collision integral is 
quite natural since there is spin fiip. On the other hand, 
the above form of T^^ implies an effective vector spin flip 
probability u'p p/ that depends on the magnetization di- 
rection n{r,t). Note that there is a priori no reason 



B. Zhang-Li relaxation time approximation 

In order to recover equations of motion for the macro- 
scopic quantities that correspond to the relaxation time 
approximation used by ZL, one needs three further as- 
sumptions: 

(i) 



i(<+ + <p,) ~ -^^S{ep, - ep). 



T^f J^(ep) 



and T^^^ are constants independent of the energy e 



P' 



(ii) 



2^ PiP p.p^ ' 



(iii) with assumption (i), at equilibrium, the collision 
integral of /y is no longer equal to zero at it should be. 
Therefore one needs to replace /y by g\i = f\\ — f^^ in the 
corresponding collision integral. 
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C. Extended relaxation time approximation 

Because of assumption (ii) the usual relaxation 
time approximation ignores some important qualitative 
physics. In fact it is possible to relax (ii) because for 
the two components /, /y it is not necessary to make any 
approximation to obtain their exact forms in both equi- 
librium and out-of-equilibrium situations. 

Indeed, by substituting the usual equilibrium form 
/°, it is easily verified that both X"^^ ^ and if^ 7^ 

0. Nevertheless J''f+ + Jjj'^" = and V^' = 0, thus 

the spin flip scattering does not modify the equilibrium 
properties of these two components. 

By contrast, as we show in the Appendix in out-of- 
equilibrium situation, by relaxing (ii) the contributions g 
and are qualitatively modified. For the perpendicular 
components, without approximation (i) it is not possi- 
ble to extract explicit closed forms for their distribution 
functions. Therefore in the following, as main assump- 
tion we assume that (i) is valid for the collision integral 
of the perpendicular components. 

In the following we only give the final expressions for 
the distributions and physical quantities in both equil- 
brum and out of equilbrum situations. The detailed steps 
of the calculations are described in the Appendix. 



D. equilibrium properties 

At equilibrium /'^ and J'^^^ are unmodified. For J^^, 
the previous equilibrium forms stay also valid but with a 
modified scattering time r ^ f . Thus at equilibrium the 
" p -isotropic" spin flip mechanism does not qualitatively 
change the physics of these two contributions. 

Essential modifications arise for the last two compo- 
nents fit,f2f At equilibrium they depend only on Cp, 
therefore the energy and angular integrations of collision 
integrals can be performed and the two components are 
easily found : 



1 



l + (277-f)^ " 

/2t(e) = -l + (2^sf)2/|| 



(36) 



For energy independent ij^^ the p integration is imme- 
diate and we obtain two perpendicular components to 
the "equilibrium" magnetization: a component mit — 
i-i-(2d''* ' )^ m|| collinear to — n x dtn and a component 
W2t = 2rf^mit collinear to —dtn. The modulus of this 
perpendicular magnetization vector is reduced by a fac- 



tor 



1 



, = compared to the collisionless situation. 

Once again, apart from this change of modulus, the spin 
flip has also induced a rotation in spin space of this trans- 
verse magnetization vector. Quantitatively, the ratio be- 
tween the two transverse components of the magnetiza- 



tion defines a new rotation angle tan 6'''^ = 277^^ = 

We have previously pointed out that, as back action on 
the d electron local magnetization, the component mn 
contributes to the effective 7 termli. The new com- 
ponent —Tsdfn2tdtn contributes to the effective Gilbert 
damping by a term 



0.21 



Tsd- 



2rj 



,sf 



mil 



(37) 



(2r,d) 



This expression exactly coincides with the ZL result that 
was calculated using the macroscopic equations of mo- 
tion in the relaxation time approximation. Taking only 
the leading order term in Tsd/T^^ we also recover the ex- 
pression of Kohno et al when the appropriate changes of 
notations are made. 



E. out-of-equilibrium properties 

To improve readability, the transport equations and 
the explicit calculation of the out equilibrium distribu- 
tions of the different components are described in the 
Appendix. The resulting physical properties are the fol- 
lowing. 

For the particle and parallel spin currents, to first order 
in p^^ we obtain : 



(38) 



Note that, even if there is no difficulty to obtain the 
results to any order in p^^, the leading order term con- 
tains all the important qualitative modifications. These 
expressions show that the spin flip scattering modifies 
both the particle current and parallel current in a dis- 
tinct manner. Such a result cannot be obtained using 
the standard relaxation time approximation and can only 
be introduced phenomenologically within the ZL macro- 
scopic approach. We remark that, when P ~ I (fully 
polarized case) , the parallel and particle current coincide 
as they should. When P = (unpolarized) the parallel 
component is zero because p^^ is implicitly proportional 
to P (see the definition above). For later use, we also 
note that to first order in p^^ the relation between the 
parallel spin current and the charge current je — ej is: 



(P+p^f-(l-P2))j, 



(39) 



As explained in the Appendix, in out-of-equilibrium sit- 
uation the spin current has also two transverse compo- 
nents collinear to dtn and n x dtn, respectively. These 
two transverse components of the spin current are easily 
accessible with our formalism but cannot be explicitly 
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found with the ZL approach and need to be calculated 
within the linear response or Keldysh approach. 

The last two components 5i,2r(p) of the out-of- 
equilibrium spin distribution function, as explained in 
the Appendix, can be recast in tensor form as gi^2i = 

T^^^Ip^P' x-i_^2{e) + (^'-^ 2/1,2 (e)], with modified functions 
2^1,2 (e) and j/i.2(e) as compared to the spin-conserving 
scattering case. With these expressions, by performing 
either angular integration or energy integration, we can 
calculate the modified energy resolved magnetic densities 
TOir(e), TO2r(e) (scc Appendix) and angular resolved mag- 
netic tensors mir(p), m2r(p) respectively, and then by 
further integration the resulting physical magnetization 
components mir,m2r. For comparison, we obtain the 
following expressions for the modified angular resolved 
magnetic tensors mir(p), m2r(p) : 

mir(p) = {^^)rn^ + <) {S'^ - 3p'p^) 
m2r(p) = {^^rn^ + Vrnf) {S'' - 3p'p^) 



with 

sf 

m I = 



Performing the angular integration wc finally obtain the 
two transverse magnetization components: 



1 



mir = 

m2r = 



1 + {2r^^^f '■'■^ 

1 + (27?«f)^ ^'"^ 



Jll. 



(41) 



As back action on the localized d electron magnetization, 
the component mir gives rise to the first spin torque term 
and thus determines the parameter u in terms of the 
charge current density je- 



sf2 



(2r,d) 



t> sf p ) eAf, ■ 



(42) 



The component m2r gives rise to the second spin torque 
term and thus determines the parameter /? as 



(43) 



This value of /? coincides with ZL result (it therefore also 
coincides with the Kohno result for spin-isotropic spin flip 
when appropriate modifications of notations are made). 
Note that our result for the value of u cannot be explicitly 
calculated with the ZL macroscopic equations of motion. 
To leading order in the ratio P/a2t — Mg/m\\ > 1 
and thus, at first sight, it seems compatible with micro- 
magnetics and experiments. We note however that other 



sources of dissipation (spin lattice relaxation for example) 
might give some contribution to the total effective Gilbert 
damping parameter a and thus modify the [3 /a ratio. It 
would therefore be interesting if experiments could pro- 
vide measurements of the different contributions to the 
Gilbert damping parameter a. 



SUMMARY, DISCUSSION AND 
PERSPECTIVES 



Using the Landau-Sillin approach, we have studied the 

transport of electrons in the presence of an effective Zee- 
man field that has a space-time varying direction. The 
key ingredient is our Ansatz form of the spin density 
matrix that consists in a linear decomposition on quasis- 
tationary distribution functions along each possible di- 
rection provided by the first order space-time gradients 
of the magnetic field direction. We have shown step by 
step how the form of the different components of the 
distribution function is affected by the presence of spin- 
conserving and spin flip scattering, in both equilibrium 
and out-of-equilibrium situations. For spin-flip scatter- 
ing we have defined a spin-basis-invariant collision inte- 
gral and an extended relaxation time approximation that 
show the mixing of the particle and parallel components. 
Our calculations also illustrate the striking difference be- 
tween a macroscopic quantity such as the transverse (per- 
pendicular) magnetization component and its underlying 
distribution. This is particularly clear for the compo- 
nents m2r and /2r(p) in the presence of spin- conserving 
scattering in the out-of-equilibrium situation. The term 
/2r(p) is non zero and leads to a finite energy resolved 
magnetization density m2r(e) when only an angular in- 
tegration p is performed. But further energy integration 
of this density gives a zero contribution as expected from 
the spin conservation rule. More qualitatively and phys- 
ically, we have clearly explained the existence of an equi- 
librium spin tranverse current in the direction n x Vj.n 
when there is no scattering at all. We have shown that, 
within the Boltzmann approach, the modulus of this 
equilibrium spin transverse current and its direction are 
affected by spin-conserving and spin flip scattering. The 
rotation is by the angle 9 and in the direction Vrn, com- 
pared to the collisionless situation. Although we have 
not calculated it explicitly, we have shown that the out- 
of-equilibrium contribution to the transverse spin current 
is purely in the direction provided by the time derivative 
n X dtn in the spin-conserving situtation. In presence 
of spin flip this contribution is further rotated in a com- 
plicated manner towards the direction dtn. Concern- 
ing the transverse magnetization components, the situa- 
tion is somewhat reversed. In "equilbrum" the transverse 
magnetization is along n x dtfi in both collisionless and 
spin-conserving scattering situations. In presence of spin 
flip it is rotated by an angle 6^^ to the direction dtn. 
For the localized d electron magnetization this rotation 
leads to the appearance of an effective Gilbert damping 
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correction term a2t- In out of equilbrum situation, for 
spin-conserving scattering the transverse magnetization 
is along n x VrU only and its modulus is directly pro- 
portional to the parallel spin current. In the presence 
of spin flip, this transverse magnetization is rotated by 
an angle 6^^ in the direction Vrn. For the localized d 
electron magnetization these two components lead to the 
two spin transfer torque terms and thus allow determin- 
ing the two parameters /3 and u. Our results for /3 and a2t 
exactly with those of the ZL macroscopic approach and 
of Kohno et al. in the spin-isotropic spin flip scattering 
situation that we consider. Our approach allows to cal- 
culate explicitly the spin current polarization Pj — u/ue 
where Ue is the electron drift velocity that is purely phe- 
nomenological in ZL macroscopic approach. 

We believe that this paper can be extended in at least 
four directions^. 

(i) For spin flip scattering, as we already pointed out, 
there is no reason to prevent terms proportional to the 
gradient of n(r, t) in the spinor probability Wp,p' (r, t). A 
natural extension of our work would consist in consider- 
ing such terms by assuming some specific forms for the 
corresponding probability Wa.p.p' {a — li, Ir, 2i, 2r). In 
fact we believe that wi,2r,p,p' could be at the origin of an 
out-of-equilibrium spin Hall current. 

(ii) In our paper we have only considered the most 

2 

simple quadratic dispersion relation gp = ^ and a p 
independant effective Zeeman field. In the spirit of the 
work of J. Zhang et al^, it would be interesting to ex- 
plore how far the equilibrium and out of equilbrum prop- 
erties are changed for a general dispersion relation and a 
p dependent effective Zeeman field. 

(iii) It is possible to adapt our method to spin-valve 
systems^. The key point consists in expanding the 
spinor distribution function /^(p, r, i) in each region 
£ — L, R,C {L = left thick magnetic layer, C = central 
non magnetic layer and R = right thin magnetic layer) 
into the most general basis to first order in time gradi- 
ents of the two magnetic layers directions nL ji{t). If for 
simplicity one assumes that the "thick layer" direction 
rii is time independent, naively our Ansatz distribution 
for the thin layer would be: 

/^(p,r,i)= /«(p,r)/ + /^(p,r,i). a, 
with 

/^(p,r,t)= f^{^,v)nR{t) 

+ (/i^(P,r)"ii X +/2^(p,r)nL) 
+T-id(/5(P. X dtriR + fi{p, r)dtnR) 
+T-i^/3?(P. X dtriR 

(44) 

In fact, preliminary results show that this Ansatz with 
an extended basis of six vectors is still insufficient^''. 

(iv) In the context of ferromagnetic Fermi liquids, the 
Landau-Sillin approach has been used for a long time 
and equations of motion of the magnetization with terms 
similar to the spin torque terms have been established for 
example by Leggett^^. Nevertheless, there are still many 
questions that concern the transverse properties. To our 



knowledge, most of the parametrizations used to study 
these systems22, are similar to that of Kohno et al. and 
Tserkovnyak et al.. We thus believe that new insights 
can be provided by a parametrization similar to ours. 

As a final remark and anticipating on our paper us- 
ing Keldysh Green functions, we have pointed at several 
places that a collision integral is intrinsically a quantum 
object and that, therefore, quantum corrections might 
affect the results. The microscopic derivation of a col- 
lision integral requires the calculation of a self energy 
which itself depends on the Green function. The quan- 
tum object that plays a role similar to the spinor dis- 
tribution f{r,p,t) is the time-space Wigner transform 
spinor Green function G(p, uj, r, t). Very similarly to the 
spinor distribution /(r, p, t) we thus propose to adopt 
the following Ansatz form for the Green function: 

G{p,uj,r,t)^ G(p,w)/ + G(p,w,r,t) • & 
with 

G(p, UJ, r, t) = Gil (p, w)n(r, t) 

+isd [Gir(p, uj)n X Vrn + G2r(p, cj)Vrn 
+Tsd [Git(p, w)n X dtTi + G2t{p,i^)dtn] . 

(45) 



VI. APPENDIX 

A. Spin-basis invariant collision integral for 
spin-flip scattering 

Spin flip scattering in a uniform ferromagnet (n(r, t) = 
no) corresponds to an "interband" process. The collision 
integral of each eigen-spin distribution function f± (p) — 
/(P) ± /||(P) is then: 

^Jdr'w^^l, [/^(epO(l-/±(6p)) 

-/±(ep)(l-/T(v))] (46) 
= Idr'nf^l, (/^(6pO-/±(6p)), 

with Wpp/ = '^|)'p'(£p')'^(^p' ~ ^p)- The corresponding 
collision integrals for the particle density and parallel 
magnetization components are: 

I^f ^Jdr' + <-,)(/(p') - /(P)) 

+^(<t-<.0(./ii(p') + /ii(p))], 

H'-fdr' -[i«^.-<';;0(/(p')-/(p)) ^ ^ 
-K<i'+<p')(/||(p') + /||(p))]- 

To extend the above results to a non uniform magnet we 
need first to find a spin-basis-invariant formulation such 
that for the uniform ferromagnet we can write a spin- 
matrix collision integral X = XI + 1- <t directly in terms 
of a spin-matrix distribution function /(p) = fl + f ■ o" 
(with / = /||no for the uniform case). The following form 
of X appears to be compatible with the above results for 



11 



T and T\\ 



imation the equations of motion become: 



J«f = /dr'<p,[/(p')(l - /(P)) - /(P)(l - Kv'))l 
with 



sf+ , sf- 

p.p' + ^'p.p 



1 I sf+ 
2(«'p,p' 



W. 



P,P 



<,P' = [sH 

and 

/(P) = f{p)I - f\\{p)no ■ d- . 

(48) 

On the one hand, the fact that /(p) appears in the colli- 
sion integral is quite natural since there is spin flip, on the 
other hand the above form of T^^ implies an effective vec- 
tor spin flip probability p, that depends on the mag- 
netization direction. Wc now phcnomcnologically extend 
the above form to space-time dependent magnetization 
with appropriate antisymmetrization to prevent imagi- 
nary terms: 



with 



Jdr' 



<p' (r. t), [f{p', r, i), 1 - /(p, r, t)]+ 



p,p' 



<p'(r,t),[/(p,r,t),l-/(p',r,i)]+ 



<,p'(r,i) = [i( 



w 



lKl-<:p'Hr,t)- ^ 



(49) 



eEVp/o = -kg - P^'^gw) + ^ / €igip') +p'''^Mp')) 

eEVp/0 



-ig\\-p^{9) 
1 1 

eEVp/i^t H -I- Wsdfif_Lti = --.g2t 



5(P')) 



+^i€(^-_^^,)Mp') 

eEVp/?, - a;«d52t = ~gu + ^ / ^(1 - 2;^)fli*(p') 
eEVp/^r + ^(^5|| - ^Vp5) +a;«d5ir = -^52r 



^Itii-^_^)g2r{p') 

eEVp/fr - u;,,52r = -^Sir + ^ / - 2;^)fllr(P') 

(51) 

The first two equations show that in this extended re- 
laxation time approximation g and (/y are both linear 
combinations of Vp/° and Vpf^. As a consequence in 
the first four equations the terms with angular integra- 
tion over p' do not contribute. To solve these equations 
it is convenient to further split each Qa into 5a = 5q + 
where ga is obtained, from the expressions found in the 
spin conserving case, by the substitution r ^ f . 



B. Transport equation and out of equilbrum 
distribution functions in the presence of spin flip 
scattering 



Using the form of the collision integrals introduced 
in the previous section, the equations of motion of all 
components in the presence of both spin-flip and spin- 
conserving impurities now read: 



eEVp/ 
eEVp/ii 



■ I- 



_Xsf~ 



eEVp/2r + J^(^/|| - ^Vp/) + U^sdflr = l2r ' 1^+ 
eEVp/lr — LOsdf 2r = ^Ir — ^^Ir"*^ , 

eEVp/2t + ^/|| + ujsd.f±ti = I21 
eEVp/it — ujsdf2t = lit 



7-sf-i- 

-^2t ' 



(50) 

with: X^f± = !dT'\{w^l%, ± <;,)(./(p') - /(P)) and 
= SdT'\{w^^%, ± <p,)(/a(p') + /a(p)) for a =11 

,1,2. 

The equilibrium properties are easily found and were 

described in the main text. For the out of equilbrum 
properties, within our extended relaxation time approx- 



To simplify the expressions and since it does not quali- 
tatively modify the results, we only retain the first order 
contribution in p^^ although there is no difficulty to find 
the exact forms. For the first four components, we ob- 
tain: 



5^f(p) - -p-f!^reEVp/o 
. 

5f (P) = -p^^^feEVp/o, 

^ - „sf 



gftip) = i+fi 

gftiP) = T^reE 



i^reE 



(52) 

Each of these components gives rise to some finite current 
only (particle or spin current). 

For the last two components 51, 2r we remind that the 
terms with angular integration over p' contribute. Once 

e\E\ 

again we write .91,2^ = 2^1,2(6) -|- (5*-'yi,2(e)] 

and split accordingly xi^2, 2/1,2 into xi^2 = 5^1,2 +^^2 ^^'^ 



2t] 



sf 
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yi,2 



2/1,2 + J/i^2- then obtain: 



y2 



with 



3 



3 1 
-(yi - 



1 



1 + (2r;-f) 



n2 



[2v'Ky2 + f) + (2r7^')'(yi + f ) 



X2_\ 

3 / 



(2r7-f)^ 
f)-(l + r72)(£l! + /f) 



f sf _ si 'I 1 1 I -cO ft^sd a fO\ 

J±- P -^T+n^ljU 2-<^^.'||i' 



(53) 



With these expressions, by performing the angular inte- 
gration we obtain the following modified expressions for 
the energy resolved transverses magnetization densities 
mir(e),TO2r(e): 



TOir(e) 



m2r(e) 
with 



-i(m||(6)+p^/2^ne(e))), 
= -m^ + 2r;^W(e), 



(54) 
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Mathematically, the extended basis constitutes a linearly 
independent basis in the space of adiabatic functions 
a(p, r, t) defined by a(p, r, t) = a\n{r, t) + a^n x Vrn + 
+a3Vrn + a4n x dtn + a^dtn, where the aa{p) are function 
of p only; i.e a(p, r,t) — implies aa(p) ~ for all a and 
a{p,r,t) =b{p,r,t) implies a^ip) =ba{p)- 
For spin conserving scattering we can generalize our work 
in two ways: (i) define distinct intraband scattering prob- 
abilities w"^ , for energies e±. (ii) define anisotropic scat- 
tering pro^jcSjilities Wp p/ . 

(i) In establishing the parametric forms /i,2 we do not need 
the explicit forms of /" and ./||'(ep). We only need that 
/"(tp) and f^^{ep) depend only on ep. But in fact quite 
generally we know that /'^(tp) = ^{A+{£p) + A-{£~)) and 
/||(^p) ~ l(^|| + (^p) +^||-(^p))- Clearly the components 
/i,2(ep) can also be split in two components j4i,2±(ep). 
Using this property, an immediate generalization of the 
above results consists in separating the collision integral 
into la = Ia~^ + 1a~ (with obvious notations) and to 
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specify distinct probabilities w^^, = _^ ±^(fp ~ ^p')- 

In that way, for isotropic scattering, all the previous ex- 
pressions are valid when interpreted for each index ±, and 
this allows considering the possibility of two distinct scat- 
tering times T±. 

(ii) Anisotropic scattering consists in writing Wp,p' — 
l^f w)fpPf(p • p') where Pe{p ■ p') are Legendre polynomi- 
als of order £. In the equilibrium situation, the only mod- 
ification of our results consists in replacing the isotropic 
scattering time 1/r = w°(e)i/(e) by the transport scat- 
tering time l/rtr = («'°(e) — ^w^ (e))iy(e). In the out-of- 
equilibrium situation the analysis becomes more subtle and 



is beyond the scope of the present work. 
F. Piechon and A. Thiaville, in preparation. 
J. Zhang, P.M. Levy, S. Zhang and V. Antropov Phys. Rev. 
Lett. 93, 256602 (2004), ibid 96, 019708 (2006) and J.C. 
Slonczewski Phys. Rev. Lett. 96, 019707 (2006). 
It appears in fact that, to describe a time dependent effect 
to linear order, it is necessary to consider an extended ba- 
sis of the nine following vectors nL,nR,nL x ur and the 
six new vectors dtriR^riL x dtriR^UR x dtnR,nL x {ul x 
dtnR),nL X {ur x dtriR) and ul x {ul x {ur x dtriR)). In 
this basis the mathematical properties described in^^ are 
valid. 



